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ABSTRACT 

We here present our first attempt to use Globular Clusters as tracers of their parent 
galaxy's formation history. Globular Cluster Systems of most early-type galaxies fea- 
ture two peaks in their optical colour distributions. Blue-peak Globular Clusters are 
generally believed to be old and metal-poor. The ages, metallicities, and the origin of 
the red-peak Globular Clusters are being debated. We here present our analysis of the 
ages and metallicities of the red peak Globular Clusters in the Virgo SO NGC 4570 
using deep Ks-himd photometry from NTT/SOFI (ESQ program ID 079.B-0511) for 
the red-peak Globular Clusters in combination with HST-ACS archival data to break 
the age-metallicity degeneracy. 

We analyze the combined 5, z, and Ks spectral energy distribution by compar- 
ison with a large grid of GALEV evolutionary synthesis models for star clusters 
with different ages and metallicities. This analysis reveals a substantial population 
of intermediate-age (1-3 Gyr) and metal-rich (« solar metallicity) Globular Clusters. 
We discuss their age and metallicity distributions together with information on the 
parent galaxy from the literature to gain insight into the formation history of this 
galaxy. 

Our results prove the power of this approach to reveal the (violent) star formation 
and chemical enrichment histories of galaxies on the basis of combined optical and 
near-infrared photometry. 

Key words: galaxies: individual: NGC 4570 ~ galaxies: star clusters - galaxies: 
lenticular - galaxies: formation. 



1 INTRODUCTION 

Globular Cluster (GC) systems are now recognized 
as powerful tracers of the i r parent galaxy's forma- 



tion history ijWest et"al] |2004| : iFritze-v. AlvenslebenI |2004 
iBrodie fc Stradeil 20061 ) . From their age and metallicity dis- 
tributions one can reconstruct the parent galaxy's (violent) 
star formation and chemical enrichment histories all the way 
from the very onset of star formation in the Early Universe 
to the present. 

Most early-type galaxies sh ow bimodal c olor distribu- 
tions for their GC systerns (e.g. jGebhardt fc Kiss ler- Patig 

m 



ll999l : lKundu fc Whitmor j2001al lb nPeng et al.ll2006l ): A uni- 
versal blue peak and a red peak for which the color and 
height relative to the blue peak vary from galaxy to galaxy. 
The blue peak GCs are generally agreed to be old and metal- 



* Based on observations made with ESO Telescopes at the La 
Silla Observatory under programme ID 079.B-0511 
t E-mail: r.kotuUaOherts. ac.uk 



poor, the properties and origin of the red peak GCs is still 
unclear. Scenarios for the formation of the red peak GCs 
range from in situ formation of a secondary more metal-rich 
population of GCs within their parent galaxy shortly after 
the first one (iForbes et al.lll997^ to major gas-rich mergers 
fe.g. lAshman fc Zepj|l99 ^) and hierarchical accretion events 
involving enoug h gas to trigger the formation of new GC 
populations (e.g. iBeaslev et aLll2002l ). The age distributions 
of the secondary GCs predicted by these different scenarios 
are different: almost as old as the old and metal-poor blue 
peak GCs but more metal-rich than those in the first case, 
of some intermediate age reflecting the time of the gas-rich 
merger in the second case, and broad or multi-peaked for 
a series of hierarchical accretion events involving gas in the 
third case. The metallicities of the secondary, and eventu- 
ally any further generations of GCs should reflect the ISM 
abundances in the merging or accreted objects at the time 
of merging or accretion. They could, at most, be somewhat 
higher for those GCs that formed late enough in a burst to 
already incorporate some enrichment during the burst itself. 
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Hierarchical accretion without gas and the formation of new 
generations of GCs cannot explain the red-peak GCs since 
dwarf galaxies known so far contain old and metal-poor, 
hence blue GCs. 

To determine the ages and metallicities of the red-peak 
GCs in one of those early-type galaxies with clear bimodal- 
ity in its optical GC colour distribution is the aim of our 
present investigation and should help constrain the forma- 
tion scenario for the red-peak GCs in this particular galaxy. 

1.1 Our approach to lift the age-metallicity 
degeneracy 

Optical data alone do not allow to disentangle ages and 
metallicities: Colour-to-metallicity transformations have to 
assume an age, while colour-to-age transformations are only 
valid for one metallicity. The degeneracy, however, can be 
broken by including near-infrared data that are more sensi- 
tive tojn«taJlicity rather than to age. 

lAnders et al.l (|2004l ) used extensive artificial star clus- 
ter tests and showed that observations in three passbands 
for GCs in dust-free E/SOs (or four passbands for young 
star clusters in dusty environments), spanning as wide as 
possible a wavelength-basis {U through K) and including at 
least one NIR-band (e.g. H or K) with photometric accu- 
racies 0.05 mag in the optical and 0.1 mag in the NIR 
allow to disentangle ages and metallicities and to determine 
individual GC metallicities to < ±0.2 dex, and ages 
to < ±0.3 dex. I.e., these data allow to distinguish 7Gyr 
old GCs from those ^ 13 Gyr old. 

Similar studies also using NIR-data to determine 
ages and metallicities of globular clust ers have only been 
done for a few galaxi e s until now (iPuzia et alT l2002l: 



Kisslcr-PatiE ct al.' 2002'; 'Hemp el etal] l2003l : iLarsen et all 
2005; Hcmpcl ct al. 2007). More than half of these galax- 
ies were found to host a population of GCs that is younger 
and/or more metal-rich than the ubiquitous old and metal- 
poor GC population. 

Some of these previous studies discussed average prop- 
erties of the blue and red peak GCs by investigating the 
mean colours of blue and red GC. We here derive ages, 
metallicities, and masses for every individual cluster with 
{g, z, K) photometry. This also enables us to study whether 
there is more than one generation of GCs in the red peak, 
to look for correlations of these parameters, and investigate 
their spatial distributions. 

We have chosen the Virgo cluster SO galaxy NGC4570 
(= VCC1692) for our pilot study because it has archival 
0— and g— band data from the ACS Virgo Cluster Project 
( Peng et al]|2006l ) that show two very clear peaks in the op- 
tical g — z colour distribution of its GC system. The red peak 
is about two thirds the height of the blue one, both have very 
similar widths. In particular, the optical colour distribution 
shows no evidence for substructure within the red peak or 
for a third peak. The mean g~ z colours of the blue and red 
peaks are 0.88 ± 0.01 and 1.38 ± 0.03, respectively, with a 
fraction of 39 % of all the 122 GCs detected belonging to the 
red peak. NGC4570 has Bt = 11.82 mag (Cote et al . 2006), 
which at a kine matic distance modulus of BDM = 31.16 
l|Mei et al.ll2007h gives it an absolute Mb = —19.3 mag, i.e. 
it is an average luminosit y SO. 

Ivan den Bosch et al.l (|l998l ) and 



Ivan den Bosch fc EmsellemI (|l998l ) detected a nuclear 
stellar disk with a radial extent of « 7 arcsec. From both 
spectroscopy based on II/3 and [MgFe] line indices as well 
as photometry in U, V and / they estimate an age of the 
stellar population inside that structure of ^ 2 Gyr and a 
metallicity close to solar. 

From the s ame ACS Virgo Cluste r survey that reported 
the GC colours iFeirarese et al.l (120061 ) show that NGC 4570 
shows a nested-disk structure composed of two morphologi- 
cally distinct inner and out disks. Detailed isophotal analysis 
reveals a blue stellar ring with radius 150 pc, but less than 7 
pc wide, around the nucleus. This ring leaves a clear imprint 
on the major axis g— band surface brightness profile (cf. Fig. 

200(1'). The structure of this inner re- 



104 of lFerrarese et all 



van den Bosch fc EmselleinI (|l998l ) in 



gion was discussed by 
the context of secular bar evolution. 



2 MODELS 

We used our GALEV evolutionary synthesi s models for star 
clusters (|Schulz et al.|[200 2: Anders & Fritz e-v. AlvenslebenI 
|2003l ) to compute a large grid of models for five different 
metallicities —1.7 ^ [Fe/H] ^ +0.4 and ages between 4Myr 
and 16 Gyr in time-steps of 4 Myr. Our m odels are based on 
Padova isochrones and a lSalpeteJ (|l955l ) initial mass func- 
tion (IMF) with a lower-mass limit of Miow = 0.10 Mq and 
upper mass limits between Mup « 50 — 70 Mq depending 
on metallicity. 

Since early-type galaxies in general, and NGC 4570 in 
particular, do not contain significant amounts of dust we 
did not include internal extinction into our grid but assume 
E(B — V) = throughout. We therefore only need three 
filters (HST F475W, F850LP and SOFI K,) to determine 
all relevant parameters (age, metallicity, and mass) for each 
cluster. 

Note that we do not depend on color transformations 
from the HST to Standard Johnson filters. For optimal ac- 
curacies, our models first compute spectra as a function of 
time, that later are convolved with the corresponding filter 
curves, in this case for the HST F475W and F850LP filters 
and the SOFI Ks filter, to yield the final magnitudes. 

Figure [1] shows examples of Spectral Energy Distribu- 
tions (SEDs) for star cluster models at two different ages of 1 
and 13 Gyr, and all five metallicities for each of them. The 
SEDs shown have been scaled to match the g-band mag- 
nitude of a solar metallicity star clusters with initial mass 
10^ Mq. Fluxes in all filter bands of course scale with clus- 
ter mass. At both ages the lines representing the different 
metallicities split up towards both longer and shorter wave- 
lengths. We exploit this to separate the effects of age and 
metallicity. 



2.1 AnalySED 

To obtain accurate ages, metallicities and masses for the 
detected clus ters we used t he SED analysis tool AnalySED 
developed bv lAnders et al.l ((2004) . This program compares 
the observed SED of a star cluster with a large grid of model 
SEDs and finds the best-fitting match to the observations. 
In this process ages and metallicities for each cluster are de- 
rived from the observed spectral energy distributions. Once 
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Figure 1. Spectral energy distributions (SEDs) for star cluster 
models at two different ages of 1 Gyr at brighter and 13 Gyr at 
fainter absolute magnitudes, respectively, and 5 different metal- 
licities ranging from 1/50 Zq to 2.5^0 in each case. All SEDs 
have been scaled to match the g-band magnitude of the solar 
metallicity model with initial mass 10® Mq. 



those parameters are found the observed brightnesses are 
translated into masses. 

Since we keep the full distribution and not only 
the best fitting solution we are able to determine the la 
uncertainties for all our derived parameters by summing 
up the normalized probabilities sorted from the highest to 
the lowest values until an integrated probability of 0.68 is 
reached. The uncertainty ranges for the different parameters 
are then given by their e xtreme value s reached within this 
Icr probability range fcf. I Anders et al. ( 2004'l for details of 
the method and e.g. lAnders et al l l|2004l ) for an earlier ap- 
plication) . 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 Near-infrared data 

We observed NGC 4570 in two subsequent nights (2nd to 
4th March 2007) using the ESO-NTT equipped with the 
SOFT near-infrared imager. This instrument consists of a 
Hawaii HgCdTe IKxlK chip and 0.288" pixels yielding a 
field-of-view of « 5 x 5 arcmin. To avoid non-linearities of 
the detector we chose an detector integration time (DIT) of 
6 seconds. 12 of these exposures were internally averaged by 
the readout electronics, resulting in an exposure time per 
frame of 72 seconds. Since NGC 4570 has an extent of only 
4' X 1.1', we could use one half of the chip for the object while 
obtaining a simultaneous sky-exposure in the other half of 
the chip, swapping sides after each exposure and applying 
small shifts to avoid contamination by bad pixels. 



nents: dome-flats corrected for their characteristic on-off 
pattern using the recipe provided on the NTT-SOFI web- 
site; illumination correction surface obtained from repeated 
observations of a standard star at different positions of the 
detector to remove the illumination gradient introduced by 
the dome flats; a "super-flat" created from the sky frames of 
both nights to remove remaining inhomogeneities and cos- 
metics from dust on the filter. Since the observations were 
carried out under non-photometric conditions we used the 
data from the instruments web-pag«0 to obtain the illumi- 
nation correction. 

After flat-fielding all the frames we estimated the back- 
ground in all frames by iteratively clipping values larger then 
3(T above the mean of the full frame. We then used six sky- 
frames that were obtained closest in time, scaled each of 
them individually to the sky-value of the object frame and 
subtracted their average from the object frame. We varied 
the number of frames to average, but found that six is the 
best compromise between signal-to-noise and artifacts intro- 
duced by the pupil-rotation of the alt-az mounted telescope. 

All the sky-subtracted frames were aligned by matching 
coordinates of several background galaxies in each frame to 
coordinates in a reference frame; we used galaxies instead of 
stars because due to the high galactic longitude there were 
too few stars to allow for proper matching. In a final step 
we stacked all but those frames with a sky-value deviating 
more than 2a from the average. 

Combining all those exposures from both nights results 
a total exposure time for our K-band image of toxp « 25 ks « 
9.5 hours. 



3.1.2 Photometric calibration 

We could not base our photometric calibration on standard 
stars, because during both nights observations were ham- 
pered by a varying degree of cloudiness. We therefore com- 
pared brightness profiles of the host galaxie s with calibrated 
data from the 2MASS survey fseefskrutskie et all (|2006h for 
a review) and own observations obtained later with the SIR- 
lUS instrument at the InfraRed Survey Facility at the South 
African Astronomical Observatory. Both results showed ex- 
cellent agreement within the error ranges of Amx ~ 0.03 
mag, that can easily be explained by minor differences in 
the filter transmission curves. However, we account for this 
uncertainty by adding the calibration error in quadrature to 
the photometric errors. Calibration of exposures from both 
2MASS and S IRIUS finally relied on standard stars from the 
iPersson et al.l ([l998) catalog, so our K-band magnitudes are 
in the VEGA magnitude system. 

This allows us to detect point sources down to rriKs ~ 21 
mag at the lOa-level, making the depth of our observations 
comparable to deep surveys obtained with the sa me config- 
uration, e.g. the K20 survey (jCimatti et akllioo^ . 



3.1.1 Reduction 

Data reduction largely followed the procedures outlined in 
the SOFI instrument handbook. The individual frames were 
corrected for the inter-row cross-talk and then corrected 
with a flat-field to remove pixel-to-pixel sensitivity varia- 
tions. These flat-flelds consisted of three different compo- 



3.2 HST data 

The HST data consisted of two datasets taken with the 
Advanced Camera for Survey (ACS) on-board the Hubble 
space telescope as part of the ACS Virgo cluster survey 

http:/ /www. Is. eso.org/lasilla/sciops/ntt/sofi/ 
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l|C6te et al.ll2004) . For both datasets with fihers F475W 
SDSS g) and F850LP (« SDSS z) we rehed on the On-the-fly 
reduction performed automaticaUy upon retrieval from the 
MAST ArchivtH and using the best reference files. We then 
performed an additional alignment step to ensure a match 
of coordinates in both frames as good as possible. 

Calibration of the HST data was done using the ap- 
propriate header entries from the fits-file s. For details on 
this p rocess see the HST Data Handbook (|Pavlovskv et al.) 
I2OO5I ). To avoid any unnecessary conversion between differ- 
ent magnitude systems, we performed photometry in both 
HST filters using ST magnitudes. 

3.3 Cluster selection 

Globular Cluster candida tes were selected from the HST 
images using SExtractor (|Bertin fc Arnout3ll996l ). A valid 
detection is characterized by at least 4 adjacent pixels with 
intensities of 3a above the local background, resulting in two 
catalogs with > 1000 objects each. We then cross- correlated 
these catalogs to remove remaining spurious detections as 
e.g. remaining cosmics. 

For the ~ 330 remaining candidates we derived intrin- 
sic s ource sizes using the ISHAPE-package from BAOLAB 
iLars en 1999). This algorithm in our case convolves a King 
profile with a fixed concentration parameter c = = 30 
(or equivalently log(-p-) « 1.5) but variable radii with the 
instr umental poi nt spread function (PSF) created by Tiny- 
Tim l|Kristll2004h and determines the best fitting radius via 
X'^-minimization. We rejected all objects appearing stellar- 
like (intrinsi c radius rj < 0.2 px ~ 0.8 pc at a distance 
D « 17 Mpc (|Tonrv et al.ll200ll '). 17 objects) or too extended 
(ri > 5 px PS 20 pc, 37 objects) to be a globular cluster. 
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Figure 2. Photometric errors as function of absolute magnitudes 
for the HST filters F475W and F850LP and the NTT filter Ks. 



mostly because they were not included within or too near to 
the edge of the K-band field-of-view, and in some cases be- 
cause they were out-shined by the bright galaxy background. 

To remove last outliers we introduced a colour selection 
criterion of (gsT - zst) ^ 1.0 and (zst - Ks,vcga) ^ 4.5, 
covering the color range of our models. We further removed 
all candidates with magnitude errors ^ 0.1 mag in g and 
z and ^ 0.3 mag in Ks (see Fig. [2] for the distributions of 
magnitude uncertainties as a function of luminosity in all 
three bands). This leaves us a final sample of 63 bona fide 
GCs. 



3.4 Photometry 

For the remaining « 280 GC candidates we obtained aper- 
ture photometry using commands from the ESO-MIDAS 
package. For the HST images we used an aperture radius 
of 15 pixels and a sky annulus from 17 to 20 pixels. To com- 
pensate for the low-intensity extended wings of the ACS 
PSF we applied an aperture correction of 0.067 (0.082) mag 
to the F475W (F850LP) magnitudes according to the tab- 
ul ated Enclosed Energy values of 0.940 (0.927) determined 
bv lSirianni et al.1 (|2005l ). Although these factors still depend 
both on the position of the GC on the detector and intrinsic 
colour of the GC, the implied uncertainties are much smaller 
then the photometric uncertainties and can therefore be ne- 
glected. 

To derive photometry for the ground-based /fs-imaging 
we transformed the coordinates from the combined HST cat- 
alog into physical coordinates in the K-image and then used 
these coordinates as center positions for the photometry. 
Using SExtractor to also obtain a K-band detection catalog 
yielded less reliable results and missed many of the sources 
found with the superior resolution of the HST. We used an 
aperture of 2" and an aperture correction derived from stel- 
lar photometry within the field-of-view of (0.2 ± 0.03) mag. 
For 117 candidates we could not derive a K-band magnitude, 

^ http://archive.stsci.edu/ 



4 RESULTS 

For three selected clusters we present the detailed distri- 
butions in Fig. [3] This Fig. demonstrates that the reason for 
the uncertainties in the derived parameters are mostly rela- 
tively to isolated secondary peaks in the probability distribu- 
tion. The cluster in the upper panel has a well-defined best- 
fit age of 1.12+°;^!^ Gyr, a metallicity of [Fe/H] = -O-Sto't 
and a derived mass of {2.64:tlil) X 10^ Mq. Its value for 
the best solution with [Fe/H] = —0.4 is a factor of 10 better 
than that for the somewhat lower metallicity [Fe/H] — —0.7 
paired with slightly higher age. The cluster in the mid- 
dle panel has an age of 16. Olio's Gyr, a metallicity of 
[Fe/H] = -1.7 and a derived mass of (AMtUr) x lO'' Mq. 
It is clearly seen that the old age low metallicity solution 
has a x'^-value about two orders of magnitude better than 
that for any other metallicity. Also the intermediate age 
solution for the same metallicity has a significantly higher 
X^. The cluster in the lower panel has an age of 5.Q2'^]^gg 
Gyr, a metallicity of [Fe/H] — —1.7 and a derived mass of 
(9.961^33) X lO'^ Mq. Note the large age uncertainty for 
this cluster, that directly infiuences the derived mass of the 
cluster. Again, however, only the lowest metallicity gives an 
acceptable fit. Note that since our grid of models does not 
sample the covered metallicity range finely enough. There- 
fore only the lowest metallicity model results in a good fit 
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Figure 3. x'^ values as function of metallicity and age for 3 dif- 
ferent clusters. 



hence we do not give errors on metallicity for those two clus- 
ters. 

These three GCs represent typical cases for a very good 
fit with a fairly well constrained age, a securely old (> 5 Gyr) 
but less well constrained age, and a poorly constrained age, 
respectively. 



4.1 Ages 

Figure |4] shows the age distribution of all the GCs from 
our final sample. At old ages (age > 10^° yr) we see the 
analogues to our old Milky Way GCs formed during an early 
phase of galaxy formation. 

We also find a significant population of younger GCs 
with ages of 1 — 3 Gyr that must have formed during some 
violent star formation event. Possible candidates for such 
events are massive bursts of star and star cluster formation 
that come along with a merger of two gas- rich spirals or the 
accretion of a gas-rich companion. 

While the light grey histogram in Fig. |4] shows the 
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Figure 4. Age distribution of our detected globular clusters. 



age distribution of all our GCs, the dark histogram shows 
metal-poor GCs only. The old peak contains GCs with 
[Fe/H]= —1.7 close to the metallicity of the halo GCs in 
our Milky Way, but also more metal-rich GCs while the 
intermediate-age peak is predominantly made up by GCs 
with metallicities higher than [Fe/H]= —1.7. 

A detailed look at the 8 GCs in the log(age) bins 8.65- 
8.9 and 9.15-9.4 shows that within their error bars they 
might as well be part of the 8.9-9.15 Gyr peak. The GCs 
in the age bins 9.4-9.9, on the other hand, could, within 
their Icr uncertainties as well be part of the truely old, i.e. 
10 — 13 Gyr population. 



^,1.1 Old and metal-rich GCs 

4 out of the 10 old clusters have SEDs best described by 
solar metallicity. This, however, does not seem very plausi- 
ble, since it would require an incredibly fast global chemical 
enrichment to solar metallicity within a timescale of only 
about a Gyr. On nearer inspection, one of these four GCs 
has a large age uncertainty ranging down to 1.25 Gyrs, so 
it might well belong to our newly detected intermediate age 
population. The three remaining GCs are located within the 
extent of the host galaxy, so that their photometry might 
still be contaminated by light from the surrounding galaxy. 
This is also supported by the poor quality of their SED-fit 
with remarkably low unnormalized integrated probabilities 
< 10-1°. 



4.2 Spatial distribution 

There is no obvious difference in the spatial distributions of 
the old and younger GCs, as seen in Fig. [S] GC numbers, 
however, are too small for meaningful statistical tests. 

We note that the bulk of the GCs we detect belong to 
the red peak of the optical GC colour distribution. Only the 
very brightest GCs from the blue optical peak are detected 
in Kb and, hence, part of our sample. 



4.3 Metallicities 

The metallicity distribution (Fig. |SJ is dominated by a large 
population of 27 out of 63 clusters with solar metallicity, 
only 21 clusters have very low metallicities of [Fe/H] = —1.7 
in the range of Galactic old GCs. 

The fact that we find only relatively few old and metal- 
poor GCs can be explained as a selection effect: Since we 
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Figure 5. Spatial distribution of our Globular Cluster sample, overlayed over our K-band image obtained with NTT/SOFI. Blue circles 
mark GCs with best-fit ages < 3 Gyr, while red circles mark GCs older than 3 Gyr. Most of the bright sources are either stars or (in 
most cases) background galaxies, that are resolved on the HST images. 
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Figure 6. Mctallicity distribution of our GCs. Light grey his- 
togram: all GCs, dark: GC with ages younger than 5 Gyr. 

only include clusters that were also detected and also have 
good photometry in we prefer intrinsically red objects. 
Old Milky Way- type globulars, however, have blue colors 
due to their low metallicity and are hence not included in 
our sample. Indeed, the bulk of all our GCs, and of the 
metal-rich ones in particular, are younger than 3 Gyr, again 
a selection effect, since 3 Gyr old GCs are brighter than 13 
Gyr old ones by 1.1 to 1.3 mag, depending on metallicity 
and filter. 



4.4 Masses 

Fig.[7]shows the mass distribution of our GCs, again for the 
full sample in light grey and for the subsample of GCs with 
ages < 3 Gyr in dark grey. The mass distribution of the full 
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Figure 7. Mass distribution of the GCs in NGC 4570. The light 
grey histogram is for all GCs from our sample, the dark one for the 
subsample of GCs younger than 3 Gyr. The solid line is a Gaus- 
sian fit to our unbinned GC mass distribution with a median of 
log((MGc[M0])) = 5.3 and (t(Mgc) = 0.6 dex, normalized to the 
number of GCs in our sample. The dotted line gives the Gaussian 
for the Milky Way GCs with log{(MGC [M0])mw) = 5.47 and 
o"mw(Mgc) = 0.50 dex, normalized to the number of GCs in our 
sample. 



GC sample clearly looks like a Gaussian with a turn-over 
mass around log((MGc[MQ])) — 5.3 and cr(MGc) = 0.6 dex, 
very similar to the turn-over of the GC mass function in the 
Milky Way which occurs at log((MGcfM0l)M w) = 5.47 with 
(Mgc) = 0.50 dex (jAshman et al.lll995l '). The fact that 
the mass distribution of the young GCs does not extend to 
the same high masses as that for the old ones cannot be a se- 
lection effect, since more massive clusters of young age would 
be easily detectable if they were there. It seems that the sec- 



Young GCs in an old SO: Formation History of NGC4 570 7 




8.5 9 9.5 

log(age [years]) 



A(g,z) < 0.10. AKs < 0.20 
A(g,z) < 0.05, AKs < 0.10 



-0.7 -0.4 
Stellar metallicity [Fe/H] 



A(g,z) < 0.10, AKs < 0.20 
A(g,z) < 0.05, AKs < 0.10 



3.5 4 4.5 5 5.5 6 6.5 7 

log(mass [solar masses]) 

Figure 8. Age, metallicity, and mass distributions for the sub- 
sample of GCs with high (light grey) and very high (dark grey) 
photometric accuracies. 



ondary event that formed the ~ 2 Gyr old GCs in NGC 4570 
did not form the GCs with the same mass function as the 
old GC population. Careful modeling of cluster destruction 
effects were required to really prove this conjecture. Dynam- 
ical friction, which is most important for massive clusters, 
might have preferentially destroyed those, in particular if the 
secondary GC population were more centrally concentrated 
than the primary one, as previously found in many, but not 
all, bimodal GC systems. In addition to that we might have 
missed massive GCs in front of the bright galaxy background 
if they either were initially closer to the centre (mass segre- 
gation) or driven towards the centre by dynamical friction 
faster than the lower mass ones. 



4.5 The subsample of GCs with high photometric 
accuracy 

In Fig. |8] we present our results for the age, metal- 
licity, and mass distributions of the subset of 
GCs from our sample with the best photomet- 
ric accuracies: Ag, Az<0.10, AKs < 0.20 mag and 
Ag,Az < 0.05, AKs < 0.10 mag, respectively. This Fig. 
shows that our results are unaffected by the slightly higher 
uncertainties in our full sample. The age, metallicity, and 
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Figure 9. Age distributions of our GC sample under the assump- 
tion of a fixed low metallicity [Fe/H] = —1.7 (upper panel) and 
[Fe/H] = 0.0 (lower panel). Bin sizes are the same as in Fig 3] 
only the two histograms for free and restricted metallicities are 
plotted next to each other to better reveal the differences. 



mass distributions are robust. High photometric accuracies 
are only achieved for the brightest clusters. Hence it is no 
surprise that the high accuracy subsample lacks some of 
the low mass GCs. 



4.6 Alternative solutions 

4-6.1 Ages for fixed metallicities [Fe/H] — —1.7 or 
[Fe/H] = 0.0 

In this subsection we explore to what extent our results 
change if we artificially restrict the allowed parameter range. 
Our aim is to see whether our basic finding of ~ 1 — 3 Gyr 
old/young GCs in NGC 4570 is affected by these assump- 
tions. 

If we limit the available parameter range to a metallic- 
ity [Fe/H] = —1.7 (see. upper panel of Fig [9] for compari- 
son) many clusters are better described by older ages. These 
are needed to counter-balance the bluer colours of the lower 
metallicity models. The older ages, in turn, also result in 
slightly higher masses. However, the principle structure of 
the age distribution and, in particular, the peak at inter- 
mediate ages of 1-3 Gyr remain. This "worst-case" scenario 
emphasizes the robustness of our cluster age determination. 

We also determined ages for a metallicity fixed to the 
solar value and again found little difference compared to the 
original results (see lower panel of Fig. [Ojl . Although clusters 
on average get slightly younger in this case, the peak of 
the distribution around 1-3 Gyr still remains, only the tail 
towards younger ages become slightly more pronounced. 

We conclude, that our finding of a substantial popula- 
tion of intermediate-age GCs with ages of order 1-3 Gyr 
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and approximately solar metallicities in NGC 4570 is very 
robust, even under extreme test changes in metallicity. 



5 DISCUSSION 

5.1 Formation scenarios for red-peak GCs 

Our detection of a substantial population of 1 - 3 Gyr old 
GCs is most naturally explained in terms of a merger or ac- 
cretion event which involved significant amounts of gas and 
triggered a strong starburst in which, together with some 
hitherto undetected population of field stars, the presently 
observed GC formed. The star clusters we find have all the 
properties of GCs, they are compact with half-light radii of 
order 3-7 pc and have masses in the range of typical GC 
masses. Those clusters that we observe with ages of 1 - 3 
Gyr certainly are only a small fraction of all the star clus- 
ters formed in such an event. They have already survived 
the most dangerous phase in their lives, the infant mortality 
and dynamical restructuring phase after the first SNe have 
expelled the gas left-over from their formation. They proba- 
bly have also survived the violent relaxation phase that re- 
structured the merging/accreted galaxies into the presently 
observed SO. The sheer number of young GCs requires sub- 
stantial gas masses involved and their high metallicities in- 
dicate that it could not have been a minor accretion event 
involving e.g. an SMC type galaxy swallowed by some ma- 
jor gas-free/poor E/SO/Sa type galaxy. Solar metallicities in 
the ISM out of which the new GCs were formed require at 
least one Sbc- or Sb-type progenitor. Sa-type galaxies have 
higher metallicities, but lack the required amounts of gas to 
fuel the starburst. 

It is intriguing that the age we find for the young 
GCs agrees very well with the stellar population age 

of ^ 2 Gyr estimated by Ivan den Bosch fc EmsellemI 

1 19981) for the nucle a r stel la r disk in the host galaxy 



('van de n Bosch et al. 19981: [van den Bosch &i EmsellemI 

[l998; S corza fc van den Bosch! 1 19981 ). Our detection of 
a substantial population of GCs does not, however, 
support the internal bar instability scenario they favor 
for the formation of this nuclear stellar disk. Nuclear 
stellar disks are in fact seen in dynam ical simulations 
of galaxy mergers/accretion ev ents (cf. iBournaud et al.l 
12004 ISpringel fc Hernauistll2005l ) , as well as in real merge r 
remnants, e.g. in Arp 214 and Arp 224 l|jog fc Chitrell2002h . 
The fact that NGC 4570 does not feature any obvious tidal 
tails can be understood as a consequence of its location in a 
high galaxy density region of the Virgo cluster, where tidal 
tails are shred down as soon as they start to develop. 

We hence suggest that the young GCs we detected and 
the nuclear stellar disk in NGC 4570 have been formed by 
the same merger/accretion event. 



6 FUTURE PROSPECTS 

HST NIC3 data would be very valuable in order to detect 
young GC closer to the centre of NGC 4570 and see whether 
the young GC system is more centrally concentrated than 
the old one or not. Of equal value would be an extension 
of the SED to shorter wavelengths. As can also be seen 
from Figure [T] for both old and young clusters the U-band is 



both sensitive to ages and metallicity and can help discrim- 
inate between both. However, imaging in further intermedi- 
ate band-passes would only be necessa ry for other galaxie s 
containing dust (for further details see I Anders et al.[20o3 ). 
Spectroscopy, albeit challenging, would allow to confirm and 
eventually further constrain the photometric ages, metallic- 
ities and masses for these young GCs, while abundance ra- 
tios, e.g. of [a/Fe] could give further clues to the progenitor 
galaxies and to the star cluster formation scenario. 

A natural next step seem to combine our results from 
this GC analysis with a detailed analysis of the stellar pop- 
ulation across the main body of the galaxy to evaluate the 
contribution (and location) of field stars from the age range 
of the young GCs. This could either be done on the basis of 
multi-band imaging in a pixel-by-pixel analy sis of the kind 
we d id for the Tadpole and Mice galaxies (|de Griis et al.l 
I2OO3I ) or analyzing the integrated spectrum with starburst 
models as we h ave done for the ~ 1 Gyr old merger rem nant 
NGC 7252 fcf. iFritze-v. Alvensleben fc Gerhardlll994al lbh. 



7 SUMMARY 

On the basis of deep NTT-SOFI Ks— band imaging in con- 
junction with archival HST ACS deep optical imaging we 
identified a substantial population of star clusters with ages 
in the range 1-3 Gyr and metallicities around solar in the 
Virgo SO galaxy NGC 4570. All these clusters are compact 
with half-light radii in the range 3-7 pc and have masses 
of order 10^ Mq. They have successfully survived their in- 
fant mortality phase and, hence, merit to be called young 
GCs. The clusters we detect in Ks make up an important 
fraction of the red peak of the bimodal G C optical colour 
distribution reported by iPeng et al.1 (120061 ) from the ACS 
Virgo Cluster Survey. 

We performed a number of test that showed our results 
to be robust. 

The ages we find for the young GCs agree well with 
stellar population ages previously determined for the nu- 
clear stellar disk in NGC 4570. We suggest that both the 
nuclear stellar disk and the young GC population have been 
formed in the same merger/accretion event. The presently 
observed GCs certainly are only a small part of the originally 
formed star cluster population, since only a small fraction 
of clusters usually survives until ages of 1 Gyr. The sheer 
number of young GCs suggests that the merger /accretion 
event must have involved substantial gas masses, the high 
metallicity of the new GCs suggests that the gas involved 
in this event must have been enriched at least to a level 
observed in present-day Sbc galaxies. 

Our analysis has shown that GC populations are valu- 
able tracers of their parent galaxy's star formation, chemical 
enrichment, and mass assembly histories. 

Clearly further analyses, in particular of the stellar pop- 
ulation across the main body of NGC 4570, are required be- 
fore we fully understand the details of the scenario that gave 
rise to the present-day Virgo SO NGC 4570 with its nuclear 
stellar disk and young GCs. 



Young GCs in an old SO: Formation History of NGC4570 9 



ACKNOWLEDGMENTS 

We thank the International Space Science Institute (ISSI) 
for their hospitality and research support and our anony- 
mous referee for very insightful comments that helped im- 
prove and clarify this paper. We are also grateful to Ha- 
gen Meyer who obtained the calibration data at the IRSF. 
RK thanks the ESO-NTT crew, Valentin Ivanov, Alessan- 
dro Edcroclitc, and Monica Castillo, for support during the 
preparation and execution of the observations. 

This publication is based on observations made with 
ESO Telescopes at the La Silla Observatory under pro- 
gramme ID 079.B-0511. This paper is also based on archival 
observations with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research 
in Astronomy (AURA), Inc, under NASA contract NAS 5- 
26555. 



REFERENCES 

Anders P., Bissantz N., Fritze-v. Alvensleben U., de Grijs 

R., 2004, MNRAS, 347, 196 
Anders P., de Grijs R., Pritze-v. Alvensleben U., Bissantz 

N., 2004, MNRAS, 347, 17 
Anders P., Pritze-v. Alvensleben U., 2003, A&A, 401, 1063 
Ashman K. M., Conti A., Zepf S. E., 1995, AJ, 110, 1164 
Ashman K. M., Zepf S. E., 1992, ApJ, 384, 50 
Beasley M. A., Baugh C. M., Forbes D. A., Sharpies R. M., 

Frcnk C. S., 2002, MNRAS, 333, 383 
Berlin E., Arnouts S., 1996, A&AS, 117, 393 
Bournaud P., Combes P., Jog C. J., 2004, A&A, 418, L27 
Brodie J. P., Strader J., 2006, ARA&A, 44, 193 
Cimatti A., Mignoli M., Daddi E., Pozzetti L., Fontana A., 
Saracco P., Poll F., Ronzini A., Zamorani C, Broadhurst 
T., Cristiani S., D'Odorico S., Giallongo E., Gilmozzi R., 
Menci N., 2002, A&A, 392, 395 
Cote P., Blakeslee J. P., Forrarcsc L., Jordan A., Mei S., 
Merritt D., Milosavljcvic M., Peng E. W., Tonry J. L., 
West M. J., 2004, ApJS, 153, 223 
Cote P., Piatek S., Ferrarese L., Jordan A., Merritt D., 
Peng E. W., Ha§egan M., Blakeslee J. P., Mei S., West 
M. J., Milosavljevic M., Tonry J. L., 2006, ApJS, 165, 57 
de Grijs R., Lee J. T., Mora Herrera M. C, Fritze-v. Al- 
vensleben U., Anders P., 2003, New Astronomy, 8, 155 
Ferrarese L., Cote P., Jordan A., Peng E. W., Blakeslee 
J. P., Piatek S., Mei S., Merritt D., Milosavljevic M., 
Tonry J. L., West M. J., 2006, ApJS, 164, 334 
Forbes D. A., Brodie J. P., Grillrnair C. J., 1997, AJ, 113, 
1652 

Fritze-v. Alvensleben U., 2004, A&A, 414, 515 
Fritze-v. Alvensleben U., Gerhard O. E., 1994a, A&A, 285, 
751 

Fritze-v. Alvensleben U., Gerhard O. E., 1994b, A&A, 285, 
775 

Gebhardt K., Kissler-Patig M., 1999, AJ, 118, 1526 
Hempel M., Hilker M., Kissler-Patig M., Puzia T. H., Min- 

niti D., Goudfrooij P., 2003, A&A, 405, 487 
Hempel M., Zepf S., Kundu A., Geisler D., Maccarone 

T. J., 2007, ApJ, 661, 768 
Jog C. J., Chitre A., 2002, A&A, 393, L89 



Kissler-Patig M., Brodie J. P., Minniti D., 2002, A&A, 391, 
441 

Krist J., 2004, Tiny Tim User's Manual. 6.3 edn 
Kundu A., Whitmore B. C, 2001a, AJ, 121, 2950 
Kundu A., Whitmore B. C, 2001b, AJ, 122, 1251 
Larsen S. S., 1999, A&AS, 139, 393 

Larsen S. S., Brodie J. P., Strader J., 2005, A&A, 443, 413 
Mei S., Blakeslee J. P., Cote P., Tonry J. L., West M. J., 

Ferrarese L., Jordan A., Peng E. W., Anthony A., Merritt 

D., 2007, ApJ, 655, 144 
Pavlovsky M., Biretta J., Bohhn R., Chiaberge M., Cox C, 

Dressel L., Fruchter A., Giavalisco M., Gilliland R., Gon- 

zaga S., Kockcmocr A., Lucas R., Mack J., Maiz-Apollaniz 

J., Mutchlcr M., , 2005, ACS Data Handbook, Version 4.0 

(Baltimore: StScI) 
Peng E. W., Jordan A., Cote P., Blakeslee J. P., Ferrarese 

L., Mei S., West M. J., Merritt D., Milosavljevic M., Tonry 

J. L., 2006, ApJ, 639, 95 
Persson S. E., Murphy D. C, Krzominski W., Roth M., 

Rieke M. J., 1998, AJ, 116, 2475 
Puzia T. H., Zepf S. E., Kissler-Patig M., Hilker M., Minniti 

D., Goudfrooij P., 2002, A&A, 391, 453 
Salpeter E. E., 1955, ApJ, 121, 161 

Schulz J., Fritze-v. Alvensleben U., MoUer C. S., Fricke 

K. J., 2002, A&A, 392, 1 
Scorza C, van den Bosch F. C, 1998, MNRAS, 300, 469 
Sirianni M., Jee M. J., Bem'tez N., Blakeslee J. P., Martel 

A. R., Mcurcr G., Clampin M., De Marchi G., Ford H. C, 

Gilliland R., Hartig G. P., lUingworth G. D., Mack J., 

McCann W. J., 2005, PASP, 117, 1049 
Skrutskie M. P., Cutri R. M., Stiening R., Weinberg M. D., 

Schneider S., Carpenter J. M., Beichman C, Capps R., 

2006, AJ, 131, 1163 
Springel V., Hernquist L., 2005, ApJ, 622, L9 
Tonry J. L., Dressier A., Blakeslee J. P., Ajhar E. A., 

Fletcher A. B., Luppino G. A., Metzger M. R., Moore 

C. B., 2001, ApJ, 546, 681 
van den Bosch F. C, Emsellem E., 1998, MNRAS, 298, 267 
van den Bosch F. C, Jaffe W., van der Marel R. P., 1998, 

MNRAS, 293, 343 
West M. J., Cote P., Marzke R. O., Jordan A., 2004, Na- 
ture, 427, 31 



